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Abstract 

A unitary coupled-channels model is presented for investigating the decays of heavy mesons and 
excited meson states into three light pseudoscalar mesons. The model accounts for the three-mesons 
final state interactions in the decay processes, as required by both the three-body and two-body 
unitarity conditions. In the absence of the Z-diagram mechanisms that are necessary consequences 
of the three-body unitarity, our decay amplitudes are reduced to a form similar to those used 
in the so-called isobar- model analysis. We apply our coupled-channels model to the three-pions 
decays of ai(1260), 7^(1670), 7T2(2100), and D° mesons, and show that the Z-diagram mechanisms 
can contribute to the calculated Dalitz plot distributions by as much as 30% in magnitudes in 
the regions where /o(600), p(770), and /2(1270) dominate the distributions. Also, by fitting to 
the same Dalitz plot distributions, we demonstrate that the decay amplitudes obtained with the 
unitary model and the isobar model can be rather different, particularly in the phase that plays 
a crucial role in extracting the Cabibbo-Kobayashi-Maskawa CP-violating phase from the data 
of B meson decays. Our results indicate that the commonly used isobar- model analysis must 
be extended to account for the final state interactions required by the three-body unitarity to 
reanalyze the three-mesons decays of heavy mesons, thereby exploring hybrid or exotic mesons, 
and signatures of physics beyond the standard model. 

PACS numbers: 13. 25.-k,14.40.Rt, 11.80. Jy 
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FIG. 1. Af*-decay amplitude. 

I. INTRODUCTION 

It has long been recognized that hadrons lying outside of the conventional constituent 
quark model must exist within the framework of Quantum Chromodynamics (QCD). These 
so-called "exotic" hadrons, speculated as tetra-quark states or hybrid states or glueballs, 
have been predicted by various calculations using the Lattice QCD, the QCD sum rule, and 
the flux-tube model, as reviewed in Ref. Thus, quite a few experimental programs have 
been developed to search for exotic mesons via the three-mesons production reactions, such 
as irN M*N nn-nN @-[|, 7^ -> M*N ttttttN, ttKKN @, Q, and NN — >■ M* — > 
mrr] Q, where the intermediate mesons M* could be exotic. To identify M*, the main task 
is to extract the partial-wave amplitudes from the final three-mesons distributions. So far, 
this has been done mainly by using the isobar model, within which two of the three mesons 
form a light flavor excited meson R (f ,p,K*, etc.) and the third meson is treated as a 
spectator in the decays of R, as illustrated in Fig. [TJa). There, the propagation of R is 
commonly described with the Breit-Wigner parametrization or with the two-body unitary 
i^-matrix parameterizations jsl, [9| constrained by the dispersion relations. In any case, the 
three-body unitarity is missing in those analyses. The noninteracting cR amplitudes, where 
c is a spectator light pseudoscalar meson, and an appropriately parametrized nonresonant 
amplitude are then summed coherently with multiplicative complex parameters which are 
adjusted to fit the Dalitz plot of the measured three-mesons distributions. 

In an isobar-model analysis of the 7r~p — > 7t~7t~tt + p and tr~p — > 7r~7r°7r°p data from the 
E852 experiment 0-0], ai(1260), a 2 (1320), 7r 2 (1670) and a4(2020) resonances were identified, 
and the exotic J PC = 1~ + meson near 1.6 GeV [7Ti(1600)] proposed from an earlier analysis 
was ruled out. The CLAS analysis [6| of 7p — > ix + n + ix~ n data identified a2(1320) and 
7r 2 (1670), but neither ai(1260) nor exotic 7r 1 (1600) at the expected levels. On the other 
hand, the COMPASS experiment [l(| claimed to have observed 7r!(1600) in the 7r~7r~7r + 
final state from a pion scattering on a lead target. As a step to understand the differences 
between the results from these analyses as well as from the previous isobar-model analyses, it 
is necessary to first examine the extent to which the isobar model is valid. This is also needed 
for developing a theoretically sound approach to analyze the three-mesons photoproduction 
data that will be obtained at JLab with the 12-GeV upgrade 0. 

The isobar model has also been commonly used to analyze the data of three-mesons 
decays of J/ip [III, Ijj], D 13-23] and B 24-35] mesons. The B and D decays have been 



analyzed with interests in the CP violation and physics beyond the standard model. Some 
B decay processes have also been analyzed using dispersion relations (36[, neglecting the 
interactions between the outgoing two-mesons subsystem and the third meson, as assumed 
in the isobar model. The strong phases arising from the final state interactions in the decay 
processes are essential in determining the weak decay amplitudes of these heavy mesons 
and searching for physics beyond the standard model. For example, BABAR (la . I20I . 2~D 
and Belle 22, 23|| extracted the Cabibbo-Kobayashi-Maskawa (CKM) CP violating phase 
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7 from the data of B T ->■ D° (or D°) K T ->■ (A^7r+7r^)ir T . They utilized the fact that 
the interference between the decay amplitude of B T — > D°K T — > (Kg7r + 7r~)K T and that of 
B T — > D°K~ — > (Kgir + ir~)K T is proportional to e Ti7 . Clearly, the accuracy of the phases 
of the partial- wave amplitudes of D°(D°) — > (KgTT + ir~), which were determined within the 
isobar model, is crucial in extracting this fundamental parameter 7 from the data. In the 
isobar model, the strong phases from the final state interactions are partly accounted for by 
using complex D — > nR couplings. However, the phases of the amplitudes generally depend 
on kinematics and have to satisfy the three-body unitarity, which is beyond what the isobar 
model can achieve. 



The above discussions strongly indicate the need for investigating the extent to which 
the isobar model is valid. Within the well-developed three-hadron scattering models, as 
reviewed in Refs. 37-39[, the isobar model is clearly a simplification since one of the mesons 
from the decay of the propagating resonance R can interact with the third meson to form 
another R. This interaction is traditionally called the Z-diagram, as illustrated in Fig. QJb). 

and iTiiN 39N41| systems that 
as illustrated in Fig. [U^c), is 
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It was well established in the studies of nN |37| , nNN 
the multiple scattering due to the Z-diagram mechanisms, 
essential to preserve the three-body unitarity for interpreting the data correctly. Only very 
limited similar attempts have been taken recently to analyze the three-mesons decay of 
heavy mesons 42|, |43 . 



The main purpose of this work is to apply the unitary approach developed in Ref. [39 
(hereafter referred to as MSL) to investigate the importance of the Z-diagram mechanisms 
in analyzing the data of three-mesons decays from heavy mesons and excited meson states. 
We will present a model that satisfies the two-body and three-body unitarity conditions. 
We start with a model Hamiltonian defined by (bare) vertex interactions f a b,R and T cR ^m* 
and two-body interactions v c / R > tCR , where a,b,c are physical light pseudoscalar mesons (ir, 
K etc.), R is a light flavor excited state decaying to two light pseudoscalar mesons, R = 
/o, p, / 2 , K*, and M* is a heavy meson decaying to three light pseudoscalar mesons. 
The vertex interactions f a b,n are determined by fitting the empirical ab — > ab scattering 
amplitudes, and are used to define the propagation of R and to calculate the one-particle- 
exchange Z-diagram amplitudes Z c > R ^ cR . The cR — > c'R' scattering amplitudes T c > R ^ cR are 
then calculated from Z c > R ^ cR by solving a set of coupled-channels equations to account for 
the three-mesons final state interactions of heavy-meson decays. In the absence of the Z- 
diagram mechanisms, our decay amplitudes are reduced to a form similar to those used 
in the isobar-model analysis. Thus we will be able to examine the effects of Z-diagram 
mechanisms in determining the Dalitz plots and the parameters of resonances that decay 
strongly into three-mesons. The model is applied to investigate the three-pions decays of 
ai(1260), vr 2 (1670), vr 2 (2100), and D° mesons. 



The organization of this paper is as follows. In Sec. [TTJ we present our model Hamiltonian 
and describe the derivation of a set of coupled-channel equations for calculating the meson- 
R scattering amplitudes from the Z-diagram mechanisms, and how these amplitudes are 
used to calculate the three-mesons final state interactions in heavy-meson decays. The 
procedures for applying our model in practical calculations are given in Sec. IIHI The results 
for calculations of the decays of ai(1260), 7r 2 (1670), 7r 2 (2100), and D° mesons are presented 
in Sec. [TV] Summary and outlook are given in Sec. [V] 
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II. FORMULATION 



Following the MSL formulation 39| of hadron reactions, we assume that the decays of 



heavy mesons into three mesons can be described by the following Hamiltonian, 

H = H + H', (1) 

where Hq is the free Hamiltonian of the considered degrees of freedom: the bare heavy mesons 
M* = ai, 7T2, D°, the bare light flavor excited mesons R = f ,p,f 2 ,..., and the physical 
ground pseudoscalar mesons denoted as a,b,c = n, K etc. The interaction Hamiltonian H' 
is defined as (In this section, the summation runs over the momentum, spin, and isospin 
spaces of the particles.), 

ff^EE^ + r^l+fi", (2) 

M* cR 

H"= £ ^, cfl + EE[/«w + 4]. (3) 

c'R',cR. R ab 

where v c ' R r ;CR denotes the cR — » c'R' transition potentials; T cRjM * (f a b,R) is the bare vertex 
describing M* — > cR (R — > ab) processes. Here we note that the term H" does not include 
any interactions with the M* states and we have neglected the interactions between the 
particles (a, b) in the decay channels of the light flavor excited meson states R. Throughout 
this paper, we will use the "right-to-left" ordering for the channel indices. (Note that 

Fm*, c r = F c r,m* and fn,ab = flb,R for the bare vertices.) 

Starting with Eq. ([[]), the reaction T-matrix is defined by the following equation, 

where E is the total scattering energy in the center-of-mass system. Since the considered 
Hamiltonian is hermitian and energy independent, it is straightforward to show that the 
S'-matrix S(E) = 1 — 2-nibiE — Hq)T(E) is unitary S^(E)S(E) = 1. This is the simplicity of 
this formulation to have a unitary reaction model. To solve Eq. (j2J), it is convenient to first 
define a scattering equation for calculating the effects only from the non-M* Hamiltonian 
H on the scattering of the ground pseudoscalar mesons c = it, K from the light flavor 
excited meson states R = fo,p,f2, •••• Namely, we will first calculate the amplitude 

T , dRlyCR {E) = {dR'\T\E)\cR) 

where the non-M* scattering operator T'(E) is defined by 

T\E) = H" + H" E _ P R+ . H", (5) 

with H = H + H" . The intermediate states in the above equation are restricted by the 
projection operator P defined by 

P = J2\cR}{cR\ + J2\ abc )( abc \- ( 6 ) 

cR abc 
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By further applying the standard projection operator method |44j,|45(, as detailed in Ref. [39 



for a ttttN Hamiltonian, one can cast Eq. dSJ) into a form for practical calculations of 
T' C , R , cR (E). By simply changing the particle labels and dropping the contributions from 
the direct 713, v 23 , v 33 interactions in the Appendix B of Ref. |39|, we can obtain the scatter- 
ing amplitudes for this investigation. The resulting cR — > c'R' amplitudes, which describe 
the multiple scattering mechanisms followed by a M* decay as illustrated in Figs. QJb) 
and[TJc), are defined by 

T c'b.\cr( e ) = Vc'r>,cr{E) + Vc'R',c'"W"{E)G c ^ RllltCllR n(E)T^ RllcR (E). (7) 

c"'R"',c"R" 

Here the driving term is 

V c 'r',cr(E) = v c i R ^ cR + Z c / R /^ cR (E), (8) 

where v c i R ^ cR is the cR — > c'R' transition potential, and the second term is the Z-diagram 
defined with the R — > ab vertex as 

Zc'R!,cr{E) = ^/R'.cc"-^ ^ ^ ^ —fc'c",R- (9) 

„ tii — iii c — ii/ c i — Hi c ii + ze 



Here c" is the exchanged meson. We have also introduced a notation E c = \jm 2 c + to 
denote the free energy operator for a particle c with mass m c and momentum p c . 
The Green function in Eq. ([7]) is defined by 

[G- l (E)] clRI , cR = 5 d>c [{E -E c - E R )5 R ,, R - E RI>R (E - E c )\ . (10) 

The self-energy of the propagation of R in Eq. (11 01) is 

ab w - h, a - bj h + le 

where B a b is a factor associated with the Bose symmetry of mesons: B ab = 1/2 if a and b 
are the identical particles or otherwise B a b = 1. 

The self-energy ffTTj) also determines the ab — > a'b' scattering amplitudes. In the center- 
of-mass system, it has the familiar form 

T a . Viab {w) = E {B alv ) l l 2 f al y, R \d~\w)] RlyR {B ab ) l l 2 f R ^ (12) 

R,R! 

with 

[d(w)] R > >R = (w - m R )5 R > )R - E R , iR (w). (13) 

We thus can determine the mass m R of bare R state and the vertex interaction f a b,R by fitting 
the empirical amplitudes of the meson-meson scatterings such as im — > tttt and irK — > ttK. 
This then allows us to predict the Z-diagram effects on T' C , R , cR through solving Eq. ([7]). 

The transition potential v c i R r, cR can be derived from phenomenological Lagrangian by 
using the method of unitary transformation It can also be taken from more funda- 

mental modelings within QCD. This is beyond the scope of this paper, and we set v c > R i tCR = 
in solving Eq. (J7J). Thus the final three- mesons scattering effects predicted in this work are 
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only the necessary consequence of meson-meson scattering under the three-body unitarity 
condition. 

The amplitude for the three- mesons decay of M*, M* — y abc, is 



T abcM ,(E) = (* { ah ) c (E)\H'\M*}, 
where the three-mesons scattering wave function is defined by 



(14) 



*£{E!)\ = (abc\ 



1 + H"- 



E-H + ie 

with (abc\ being the three-mesons plane- wave state. From Eqs. (J2J) and (j3J), we see that 

h'\m*) = J2\cR)(cR\r cR ,M*\M*), 

cR 



(15) 



cyclic 



(abc\H"= ]T J2( abc \f^',RWR)(c'R\ 

(a'b'c') R 



cyclic 



Here the symbol V] means taking summation over the cyclic permutation, (a'b'c') = 

(a'b'c') 

(abc), (cab), (bca). Because of the orthogonality conditions, (cR\M*) = and (abc\cR) = 0, 
the above relations allow us to write Eq. (114j) as 



T abcM »(E) = (abc\H" 



P 



1 



D ' 



P 



D 



H'\M*), 



E-H + ie 

where Pd is the projection operator for the space spanned by cR states, 

P D = Y,\cR)(cR\. 

cR 



(16) 



(17) 



Following the procedures in the Appendix B of Ref . (39| , one can show that 
1 



Pi 



D 



E-H + ie 



Pd= J2 \c'R')G cmi>cR (E)(cR\ 



c'R'cR 



c'R',c"'R"' (E)T' clllR „, ,, fiii (E)Gc"R", c r(E) (cR\ , 

cR,c'R' c"'R"',c"R" 

(18) 

where T^„ R ,„ C „ R „(E) and G c r, c 'R'(E) have been defined in Eqs. (J7J) and ( |T0l) . respectively. 

Substituting Eq. ffl8l) into Eq. ffl6|) and using the vertex functions of H' defined by Eq. (j2J), 
we can write 



cyclic 

T a bc,M*(E) = ]P T(a'b')c',AI*(E), 
(a'b'c') 



(19) 



with T(afy Ci M*{E) being the amplitude for the subsequent decay of M* — y Rc — y (ab)c 
expressed as 



rp l rp\ T-iIsobar / tt<\ < t-iFSI / rp\ 

J-(ab)c,M*y^) — 1 (ab)c,M*\^) + 1 (ab)c,M*\^ ') ) 



(20) 
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FIG. 2. Amplitude for the three-mesons decays of M* in the unitary coupled-channels model. The 
bulb labeled T' is the T-matrix element for the cR — > c'R 1 process without M* excitation [Eq. J7|)]. 
The dressed cR Green function, yielding G cR:CR >(E), is indicated by the gray circle [Eq. (fTQ|)]. 



where FSI stands for final state interaction, and 



ilsobar 



= E Y,( ab \fab,RGcR,c'RiE)r c , R ,, M *\M*), 



(21) 



(ab)c,M 



R c'R' 



T(ab) c ,M*{E) ~ EE E { a ^fab,RGcR,c'R'{E)T' c , RlclllRl u(E)G c H IRII ^ c H R H(E)T c H RIIM *\M 
R c'R' c"'R'",c"R" 



(22) 



Equation ( 120]) is illustrated in Fig. [2j We now note that the commonly used isobar- 
model analysis corresponds to keeping only the term T^ 3 ^»(E) within our formulation. 
The difference between different isobar-model analyses is in the parametrization of the Green 
function G cRjC > R i(E) and the vertex function f a b,R- We also note that even in this simplified 
case, G cR ^ c i R i(E) and f a b,R are related through Eqs. (ITDjl and (TTTj) within our formulation, 
but are often not treated consistently in the isobar-model analysis. 

The full decay amplitude (120]) can be concisely written as 



Obviously, Tj$fcfa.(E) of Eq. Q2I) can be obtained from Eq. ^ by replacing \T dR , jM *) 

with r c ,R,, M *\M*). 

For strong decays, the resonance pole positions and decay widths of heavy mesons (M*) 
can be shifted by three-mesons scattering. This can be seen by considering the M* propa- 
gator defined by 



T(ab)c,M*{E) = ^2 ^2{ a b\fab,RG cR , c > R '(E)\r c > R r jM *), 
R c'R' 

where the dressed M* — > cR vertex function is defined by 



(23) 




(24) 



G M *(E) = (M*\— -L— |M*>. 




(25) 



With the projection operator method, as applied in Ref. 39[, 



one can show that Eq. (1251) in 



the rest frame of M* can be written as 



G^\(E) = E-M Mt -E M ,(E), 



(26) 



where M^-* is a bare mass and 



T, M *(E) — 2J (M*\r M * tCR G cRnC ' R /(E)\r c r R i tM *). 

cR,c'R' 



(27) 



Here \T c 'R',m*) is defined in Eq. (121 

The resonance pole positions, E po \ e , are defined as zeros of G]^,(E). They are on the un- 
physical sheets of the complex-energy Riemann surface and are thus defined by the following 
equation: 

Ga/.^poIc) = 0. (28) 



We use the analytic continuation method of Refs. [471 . |48j to solve Eq. (|28|) and find E po \ c 
for the considered coupled-channels model with unstable meson- R channels. 

If we replace the dressed vertex \r c R,M*) by the bare vertex T c r } m*\M*) in calculating 
Sm^-^poIc) of Eq. (1271) . then the solution of Eq. (1281) is the pole position of the isobar model, 
which does not include the three-mesons final state interactions. 

We now note that the Green function (126|) can be related to the excitation of M* in the 
cR — > c'R' transition amplitude T™ RI cR (E). The matrix element of Eq. (jlj) between (c'R'\ 
and \cR) states is 

T c >r< :C r(E) = T' c , R , cR {E) + T™ R , cR (E), 

where the first term has been defined in Eq. ([7]). The second term is the "resonant" part, 
and is shown to be (using the projection operator methods) 

T res ( P \ _ ( C ' R> I Fc>R> ,M* )(£m* ,cR \ cR) . . 

1,r>, c rW- E-M^-Zm^E) ' 



where 



M*,c"R" 



c"R" 



$c"R",cR + E ^c"R" c"'R"'Gc"'R"',cr(E) 



c"'R"' 



(30) 



The self-energy Hm*{E) in Eq. (|29|) is defined in Eq. ( |27|) . and can be related to the decay 
amplitude [Eq. (fT4l) ] as follows. For a derivation, consult Appendix |Aj Consider the "decay" 
width T^ll (E) of the bare M* (not the physical resonance state) defined by 

rf t *(£) = 2ttY,KE ~E a -E h - E c )\T abCtM ,(E)\ 2 . (31) 

abc 

By using the unitarity relation for T' in Eq. (jSJ), we can actually show that the right hand 
side of the above equation is 

27rJ>(£ ~E a -E b - E c )\T abcM 4E)\ 2 = -2Im[£ M , (£)]. (32) 

abc 

Equation (132]) is used to check the accuracy of our numerical calculations of the Dalitz plots 
that are calculated from T abCt M* (E) using the formula detailed in Appendix [Bl 



III. FORMULA FOR NUMERICAL CALCULATIONS 

For numerical calculations of decays of M* into three light pseudoscalar mesons, M* — >■ 
abc (Fig. [2]), we perform partial- wave expansions of the equations presented in Sec. [Ill in the 
M* rest frame. The kinematics of this decay is specified by the following: 

M*(0, S Z M „ T Z M *) R(p R , 4, t z R ) + c(p c , 0, t z c ) 

a(p a , 0, f a ) + b(p b , 0, t z b ) + c(p c , 0, tl), (33) 
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where the variables in the parenthesis for each particle are its momentum, z components of 
the spin and isospin, respectively. 

Our task in this section is to relate the partial-wave forms of all of the equations presented 
in Sec. [TT] to the basic input that will be determined by using the available empirical meson- 
meson scattering (-kit — > hit, ttK — > ttK, etc) amplitudes. In this way, the final three-meson 
scattering effects can be predicted for investigating heavy-meson decays. 



A. R — > ab decays 



In a decay R — > ab, the spin (s R ) [isospin (t R )] of the parent R state is the same as the 
relative orbital angular momentum L ab [total isospin I ab ] of the two-body ab system. Thus, 
the partial-wave expansion of the vertex function f a b.n in the rest frame of R is 

fabAQ) = (t a t:t b tl\t R t z R )Y SR ^ R (q) f^ b (q). (34) 
Here t a is the isospin of meson a and t z a is its ^-component; (jimij^TT^I JM) is the Clebsch- 
Gordan coefficient; <f is the relative momentum between a and b; f abR ab {q) is a scalar function 

satisfying f^ b ab R ab (q) = for L ab ^ s R and/or I ab ^ t R . 
We use the parametrization 

1+CW2) 



SR,L ab Vt R ,I a 



9ab,R 



1 



(35) 



.1 + (q/Cab,R) 2 . 

The parameters g a b,R and c abtR and the bare mass m R of R are adjusted to fit the empirical 
partial-wave amplitudes. The number of bare R states included in the model depends on a 
partial wave considered and the energy region covered in the fit. 



B. The 7T7T model 

We give an expression for the amplitudes of the scattering of two light pseudoscalar 
mesons in the partial-wave basis. Here, we limit ourselves to only tttt scattering because we 
consider only three-pion heavy-meson decays in this work. To fit tttt data up to invariant 
mass W = 2 GeV, we include tttt and KK channels. Then Eq. (\V2\i with total angular 
momentum L and total isospin / in each partial-wave is of the following analytic form [note 
that/^(g) = / Q L ^(g)]: 

t^M, = E flU>(v>kU E )fkU<i), (36) 

R',R 

with 

[(r LI r\E)] R/R = (E — m R )5 R ,, R - X L R [ R (E), (37) 
where m R is the bare mass of R and 

ab=«n,KK Jo E - E a (q) - E b (q) + te 

and 

f L b ab ' lab (q) = { fe 7 " 6 ^) ( if a and b are identical Particles), 

^ fab: b R ab (<l) (otherwise). 
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C. Coupled-channels equations for cR — > c'R' scattering 

For given total angular momentum J, parity P, and total isospin T, the partial-wave 
form of Eq. (J7|) for the cR — > c'R' scattering can be written as 

roo 

{c"'R"') l m,{c"R") l , l 

Here (ci?)/ denotes the ci? state with the relative angular momentum / allowed for given 
JPT; p (p') is the magnitude of the incoming (outgoing) relative momentum of the cR (c'R') 
state. The Green function can be written as 

-E)]( c "'R"') ; ///,(c"i?");» = Si>'> > i"5 c »' !C n{[E — E c »(q) — E R »(q)]6 R "r !R " 

-^ c r'",r" (?) E - -E c »(g))}, 

(41) 

where the self-energy Y< R , R (q, w) is calculated from Eq. ( ITT]) by inserting the partial-wave 

expansion (I34p and performing a Lorentz transformation to boost the function fa bR ab (q) 
from the rest frame of R to the center-of-mass frame of the cR system. We also need to 
symmetrize the intermediate states with identical mesons. Explicitly, we have 

m Rl m R r 2a M ab (q) /^S° b (?)/fe Joi (?) 



I en V / m R' m R f°° 2j 



E R ,(p)E R (p) Jo * ^[Ml b { q )+p^E-E c {p)-[Ml b { q )+pY' 2 + ie 

(42) 

where the summation is over all two-mesons states ab of R ^ ab decay, M ab (q) = E a (q) + 
E b (q). The partial- wave matrix elements Z^J R ,^ ; ( c ^) i (p',p; -E 1 ) in Eq. (140]) of the Z-diagram 
mechanisms, defined by Eq. ([9]), are given in Appendix [C] for the case that i? decays into 
two pseudoscalar mesons. 



D. The M* — > abc decay amplitudes 

The amplitude [Eq. ( ]23]) ] for a strong three-mesons decay is given by 

T( ab ) c ,M*(PaiPb,Pc] E) = J2 fab,R>(PaiPb)G cR/ , cR (p c , E)T cR , M *{p c , E), (43) 



R'R 



where the Green function G cR ' tCR has been defined by Eqs. (JJJJ) and (1421) . 

The R — 7- a& vertex function in Eq. ( 143]) is obtained from boosting the matrix element 
fahR ab {Q)> defined by Eq. (134]) in the rest frame of R, to a moving frame where p R = p a + p b '- 



fab,R{Pa,Pb) 



m R E a (q)E b (q) 
\ E R (p R )E a (p a )E b (p b ) 



(yM|t^)F Sj! , 4 (g)/^( g ) 



(44) 



where q is the relative momentum between a and b in their center-of-mass system; the 
relation among q, p a and p b can be seen in Appendix [C] [Eqs. (IC5I) and (1C6|) ]. For the strong 
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decays, T cRj m* in Eq. (143]) for the dressed M* — > cR vertex can be written down by using 
Eq. ([24]) as ' 

F C R,M*(Pc, E ) = E (^^I^M'Sa/.X^R^ 
1,1* ,si 

(45) 

where Sm* and Tap are the spin and isospin of M*, and (denoting the parity of M* as Pa/*) 

/>00 

F {cR)l , M * (Pc, E) = F (cR)l iM * (p c ) + J2 d( l V 2T {cRU™' be, q; E) 

(c"R") l „,(c'R') ll J0 

x G {c n RII)illAdR , )v (q, E)F^ B f) l , ,m« (q) . (46) 

We parametrize the bare vertex function F( cR ) u m*(p) as 

/ A 2 \ 2 +(i/2) 

W ( r\\ - 1 1 A (cR)l^* ) fP\ (A7\ 

^«,„m. w - j^^T [p> + ^U, m . J W ' (47) 

where C( cR ) h M*, ^( c r)i,m* anc ^ m T are the coupling, cutoff and the pion mass, respectively; 
A is a scale factor, and is set to be A = 1 GeV. The couplings C^ cR ) u m* are nonzero only 
when the transition M* — > (cR)i is allowed by symmetries, e.g., those are nonzero only when 
I satisfies \Sm* — s R \ < I < Sm* + sr and Pm* = Pr x (— ) /+1 . Here it is noted that for a 
strong decay the bare vertex function F( cR ) ljM *(p) is related with T cRjM * as 

r c fl,M*(p)= E (*cW£|T^t£ + 0<"*^ (48) 

For describing the weak decays of M* such as D°, the above expressions of the M* — > cR 
vertex function need to be modified to include the isospin nonconserving AT ^ transition. 
This will not be considered here. Instead, we are interested only in the importance of 
three-meson scattering after the weak decay of D°, and it is sufficient to use the above 
parametrization by extending F^ cR ) uM * in Eq. (147)) to depend on JPT of cR state. 

The amplitude for the commonly used isobar model, as defined by Eq. fl2T]) . is 

T($ c ^(p a ,p b ,p c ;E) = Y,fab,RiPa,p b )G cR ,MPc,E)T cR MPc)- (49) 

R'R 

We see from Eqs. (143]) and (149]) that three-mesons decay amplitudes of these two models 
differ from each other only in the functions describing the M* — > cR decay. The formulae for 
calculating the Dalitz plots of final three-meson distributions from these decay amplitudes 
[Eqs. (T4"3"]) and (T451) ] are given in Appendix [B] 

E. Determinations of resonance positions 

With the partial- wave expansion (145]) . the resonance pole condition G M \(E po i e ) = of 
Eq. ([28]) leads to 

E po i e = M M , + Ejvf * (-Bpole) , (50) 



11 



where [Note that F M * tic/Rf)i ,(q) = F£ c , RI)iiM ,(q)] 

E M *(E)= Y, l^dqq 2F M*,(cm>) l Xq) G (cm') l ,,(cn)Xqi E ) F (cR)i,M*(<l' E )- ( 51 ) 

(c'R') l ,,(cR) l 



c 



Here F( cR } lyM *(q, E) has been defined by Eq. (1461) ; / dq means the momentum integral is 

J c 

performed along the complex momentum path C. 

We apply the analytic continuation method developed in Refs. 47], |48| to find resonance 
poles from solving Eqs. ( |50l) and ( 15"TT) for the considered model. 



F. Relations of the isobar models with the Breit-Wigner parametrization 

For the isobar model defined within our formulation, we can establish some relations with 
the commonly used Breit-Wigner parametrization. If we neglect the Z-diagram effects by 
setting F^cR) l ,m* to F^ cR ^ l>M * and use the partial- wave expansion (1451) . Eq. (1291) for cR — > c'R' 
scattering can be written as (omitting the momentum variables) 



1 (c'R') ll ,(cR) l K^) ->• - — Q (o) ; ' 



nves,JPT fJp ^ x F (c' R') v ,M* F M* ,(cfl) 

E-M^-^l(E) 
where J = Sm* , P = Pm* , T = T M * for strong decays, and 

e£1(£)= ]T I q 2d q F M\(c>R%(q)G {c >R>) l , jiC R )l (q,E)F(cR) U M*(q)- (53) 

(c'R') l ,,(cR)i J0 

Equation ( )52|) is similar to that of the commonly used Breit-Wigner parametrization in the 
analysis using the isobar model or the i^-matrix model: 



\j & (c'R') , AnBW /ol \ P M( C R), /rBW /n] 

V«'V,( C «)A^ ~ £ — M r BW + z(r t B W/2) ' 1 J 

where r^m f is a partial decay width for M* — > (cR)i, which is related to the total decay 
width as 

r t T=ErfI )r (55) 

(cR)i 



Now we introduce s R and t R that specify the spin and isospin of R. Then Eq. (1531) can be 
written as 

E [SS^)]. r , (56) 
SR,tR 

where 

«(£)] s _ E E 

{{c! R') l ,\s R i=s R ,t R ,=t R \ {{cR)i\s R =s R ,t R =t R } 



X 



roo 

q 2 dqF M *,(c>R>) l ,{ ( l) G (c'Ri) l ,,(cR) l {qi E)F {cR)hM *(q)- (57) 
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We denote a conditional sum of (cR)i by J2{( C R),\s R =s R ,t R =t R }i i n which the (cR)i state is 
summed, keeping sr and tn constant, i.e., sr = sr and Ir = tR. Thus it is reasonable to 
make the following interpretations 



M 



BW 



■nBW 
1 tot 

SR,t R 



BW\ 



(M r BW ) 



-21m 
-21m 



Re 

,(o) 



BW\ 



BW\ 



P£l(M* w )k, 



BW\ 



(58) 
(59) 
(60) 



Equations ( l58j) -( l60j) will be used in our later comparisons with the data listed by PDG. We 
note here that the above identifications are very qualitative. 



IV. APPLICATION 



In this section, we apply our model explained in the previous sections to investigate the 
three-pions decays of heavy mesons ai(1260), 7r 2 (1670), 7r 2 (2100), and also D°. Our first task 
is to determine the parameters of our model. To simplify the calculations, we determine the 
vertex interactions R — > ab for ab = tttt, KK by fitting only the tttt scattering phase shifts 
up to the invariant mass W = 2000 MeV. This is clearly a simplification since the data 
associated with KK channel should in principle be included in our fits and we must also 
include four-pions channels that have been considered to be important in the isoscalar-scalar 
(L = I = 0) partial wave. However, such a detailed study of meson-meson scattering can 
only be done rigorously by extending our formulation to account for the direct meson-meson 
interactions v n >w,n b which must be carefully derived from effective field theory approaches, 
e.g., Refs. mm 

, to make sure that the predicted tttt amplitudes near threshold have the 
analytic properties constrained by the chiral symmetry. Furthermore, the inclusion of ?v&',a& 
in our model Hamiltonian H of Eq. will greatly complicate the scattering formulation, 
as can be seen in the ttttN formulation presented in Ref. (39[. For our present limited 
purpose of demonstrating the importance of three-body unitarity, our simplified model that 
reproduces irir phase shifts in s, p and d waves up to W = 2000 MeV is sufficient. For the 
same reasons, we neither include ttKK Z-diagrams nor make an attempt to estimate the 
errors of the determined parameters. 

In Sec. IIV A| we determine the model parameters by fitting the tttt phase shifts and reso- 
nance parameters listed in PDG. With the parameters obtained from the fit, we determine 
the pole positions (Sec. IIV B[) and calculate the Dalitz plots (Sec. IIV C|) from the M* — > hutt 
amplitudes including the Z-diagram [Eq. (143]) ] or without the Z-diagram [Eq. ( 149]) ] with the 
formula given in Appendix For the calculation without the Z-diagram, we either simply 
turn off the Z-diagram in the full calculation, or fit the isobar model to the Dalitz plot from 
the full calculation. Our main focus is to examine the effect of the Z-diagram (and thus the 
three-body unitarity) on these quantities by detailed comparison of the results calculated 
with and without Z-diagram mechanisms, thereby providing information about the extent 
to which the commonly used isobar-model analysis is valid for extracting the properties of 
heavy meson from three-mesons decay data. 
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TABLE I. Masses (M Ri ), couplings (g^^, 9kk an d cutoffs (c^^n., c K ^ R .) of the i-th bare R 
states, Ri, in the tttt partial wave with the angular momentum L and the isospin /. The couplings 
and cutoffs are defined in Eq. (j35|) . 



R (L,I) 


M Rl 






9kk,R! 


C KK,R X 


Mr 2 


9ttit,R2 




9kk,r 2 


C KK,R, 2 




(MeV) 




(MeV) 




(MeV) 


(MeV) 




(MeV) 




(MeV) 


/o (0, 0) 


1220 


-0.898 


441 


0.006 


1970 


2400 


0.700 


955 


-1.179 


394 


P (1, 1) 


891 


-0.291 


394 


0.106 


467 


1840 


0.021 


1973 


0.167 


394 


h (2, 0) 


1607 


-0.051 


567 


0.015 


818 














FIG. 3. (Color online) Phase shifts (upper) and inelasticities of the tttt scattering (lower): (left 
panels) L = I = 0, (center panels) L = I = 1, and (right panels) L = 2,1 = 0. Data are taken 



from Ref. 



5J-155]- 



A. Determinations of model parameters 



1. Fits to tttt amplitudes 

Our first task is to determine the R — >■ ab vertex function fabR ab { ( i)i defined by Eq. 
by fitting the irn phase shifts. We include hit and KK channels and use the formulae ( 155]) 
and (136|) to fit the available irn amplitude in s, p and d partial waves. In our fits, the number 
of bare R included in each partial wave is 2, 2, 1 for (sR,t R ) = (L ab , I ab ) = (0, 0), (1, 1), (2, 0), 
respectively. The resulting parameters are listed in Table HI 

As shown in Fig. [5J we are able to get good fits to the empirical partial-wave ampli- 



tudes [53H55J. The nonzero values of the inelasticities are due to the couplings to KK 
channels. For s and p waves, (L ab , I ab ) = (0, 0), (1, 1), the high quality fits are obtained only 
when two bare R states are included. It is noted that a partial-wave analysis using more 



recent data [56] has found a unique solution for W ~ 1000-1800 MeV. Although our present 
model is reasonable enough to address the question on the importance of the Z-graphs, those 
data should be considered for a more quantitative application of our model. 
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TABLE II. Pole positions of the tttt partial-wave amplitudes with the angular momentum L and 
the isospin / in the complex-energy plane. We list only the poles below Re[i£] < 2 GeV. Roman 
numerals in the square brackets specify the Riemann sheet on which the pole exists. We use the 



convention for specifying each Riemann sheet, I - IV, which is defined in, e.g., Ref. 57]. 



L 


I 






Pole positions (GeV) [Riemann sheet] 








0.43 


- 0.27i 


[II] 1.00 - 0.009? [II] 1.35 -0.17i [III] 


1 


1 


0.77 


- 0.081? 


[II] 1.61 - 0.12i [III] 


2 





1.25 


-O.lOi 


[III] 



We have also determined the resonance pole positions by applying the analytic continu- 
ation method of Refs. 47J, |48j. The results for Re [E] < 2 GeV are listed in Table [UJ It is 
interesting to note that we have two bare R states in s wave, but we have identified three 
resonance poles on different sheets of Riemann surface: sheet II is (up) consisting of the 
unphysical (u) txti and physical (p) KK sheet, and sheet III is (mm). 

We find that the poles listed in Table [TT] can be identified with the ttti resonances listed 
by PDG j58[. For the (L, J) = (0, 0) s-wave partial wave, our results can be identified with 
/o(600) (or cr), / (980) and / (1370). For (L, I) = (1, 1) p wave, our results correspond to 
the p(770) and a higher mass p. The resonance ^(1270) can be identified with our result 
for the (L, I) = (0, 2) <i-wave partial wave. Here we note that the imaginary part of the 
position (1.00 — 0.009z) GeV, which corresponds to having 18 MeV of the full width, in the 
isoscalar-scalar L = I = partial wave (the second row of Table [TT]) is too small compared 
with the full width 40 — 100 MeV of /o(980) listed by PDG. This perhaps can be improved 
only by extending our model to include four-pions channel and direct interactions v a ^ a /y 
with ab, a'b' = tttt, KK. But this is beyond the scope of this investigation, as discussed in 
the beginning of Sec. IIVI Here we mention that the imaginary parts of the pole position 
of /o(980) from some previous irir models are also smaller than the full width, (40-100) 
MeV, listed by PDG, such as 28 MeV from Ref. 0, and 29 MeV from Ref. Q. A quark 
model [6l[ also gave only 15 MeV for the width of /o(980). 

In most of the previous studies of heavy-meson decays, only the s-wave resonances are 
included as resonance poles while the p-wave poles are included in Ref. j^j. The other 
resonances are included using the Breit-Wigner form. In our calculations, we include all 
resonance poles in irir s, p and d partial waves listed in Table HTJ 

We evaluate the irR Green functions [Eq. (JJJJ)] and the matrix elements of Z-diagrams 
[Eq. fIClOj) ] with the parameters listed in Table HJ We solve the coupled-channels equa- 
tions [Eq. P0|) ] to obtain the ttR — > n'R' scattering amplitude for given JPT, Tjy^ (p',p, 
including all allowed relative orbital angular momentum between 7r and R. The resulting 
T{Sr') , 0kR)i(p'jPj-^') are then used to calculate the M* — > ixixix decay amplitudes (|4"5|) and 
find resonance poles associated with M* by solving Eqs. ( !50|) and ( 15T|) . 



2. Parameters for the decays of M* states 

To calculate the decay amplitudes for ai(1260), 7r 2 (1670), 7r 2 (2100) and D°, we now need 
to determine their bare masses Mj^*, and the parameters C^r^m*, ^-{■nR) l ,M* °f Eq. f|4T|) for 
the M* — > tiR vertex functions. Ideally, we should determine these parameters by fitting the 
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TABLE III. Properties of M* = ai(1260), vr 2 (1670), vr 2 (2100) to which our model is fitted: Isospin 
(I), spin (J), parity (P), and charge conjugation parity (C); pole masses ; branching ratios (BR). 

oi(1260) 7r 2 (1670) 7r 2 (2100) 

I(J PC ) 1(1++) 1(2"+) 1(2^+) 

Pole masses (MeV) 1230 - 213* 1672 - 130i 2090 - 313? 

BR(M* -»• vr/o) (%) 23 13 45 

BR(M* -»• vrp) (%) 74 31 19 

BR(M* -»• tt/ 2 ) (%) 2.5 56 35 



Dalitz plots of 7T7r7r distributions measured experimentally. However, such a rather complex 
process is not needed for our limited purpose here to mainly investigate the extent to which 
the commonly used isobar-model analysis is valid. It is sufficient to choose our parameters 
guided by the resonance positions and branching ratios listed by PDG [58 E The data for 
the resonances considered in this work are listed in Table IHIl 

We first notice that the data in Table II III are the averaged values from several analyses. 
Most of these analyses parametrized the M* decay amplitudes with the Breit-Wigner form, 
and all of them treated the final three pions as the paired two pions (whose correlations 
are described by either the Breit-Wigner form or the if-matrix) and the noninteracting 
spectator. We thus assume that the Tf^ /2 is the imaginary part of the pole masses in 
Table from which we can use the listed branching ratios BR(M* — > cR) to calculate 
= BR(M* cR) x If wj| The resulting values of as well as the pole masses 
of M* are then used as data to determine the parameters of F( cR ) uM *{p) [Eq. (147|) ] and 
Mjj^* [Eq. (126]) ] by using Eqs. (1501 . (I59I) and (|60|) . Because we have more parameters than 
the number of data, we use for simplicity the same cutoff for all of -F( c r) (i m* (p) for a given 
M*, and we only adjust the coupling constant C( c _R), min ,M* with the lowest allowed angular 
momentum Z min for cR; the other C^ c ^ u m* are set to zero. The resulting parameters for 
ai(1260), tt 2 (1670), and vr 2 (2100) are listed in Tables EM 

The D° meson (1865 MeV, J p = _ ) mainly decays weakly and thus three-meson scat- 
tering effects have very little effect on its mass and width. We thus will only investigate 
the Dalitz plot for the D° — > 7r+7r~7r° decay. The BABAR Collaboration [l9| presented the 
Dalitz plot data for this process, and we utilize their observation that the D° — > 7r+7r~7r° 
decay is dominated by T = np channel for simplicity. Thus we use the following param- 
eters: C( nIi ii) ltM * = 1 (see Table HVl for the notation), and C^^^m* = for other partial 
waves; A c ^m* — 1 GeV. We are interested only in the difference between the Dalitz plots 
calculated with and without Z-diagram, so this simple choice of parameters is sufficient. It 
turns out that this simple choice of the parameters well reproduces the shape of the Dalitz 
plot presented by the BABAR Collaboration jlQj ] . Clearly, the above procedure is just for 
a very rough estimate of bare M* parameters. In the future, we should fit the Dalitz plot 
data directly. But the present procedure is sufficient for our purpose in this paper. 



1 PDG lists small but non-zero branching ratios of decay channels that we do not consider in our model; 
we ignore them. 

2 In data analyses with the Breit-Wigner parametrization of the M* decay amplitudes, the partial width 
and the imaginary part of the M* pole masses are not related by = BR(M* — > cR) x rj^\ We use 
this relation just for determining the parameters with this rough estimate. 
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TABLE IV. Masses (M^*), cutoffs (AkR,m*)> an d couplings {C^ R Li^ uMt ) of the bare M* = 
ai(1260). The cutoffs and couplings are defined in Eq. (|4"7|) . For Cr nR Li\ uM *, Rf 1 means the 
i-th bare R state with the spin L and the isospin I, and I denotes the orbital angular momentum 
between R\ l and it. The second (third) column shows the parameters for the unitary (isobar-fit) 
model. The hyphens ( — ) indicate the unused parameters. See the text for the definition of the 
isobar-fit model. 

oi(1260) 



Unitary model Isobar-fit model 



M° M * (MeV) 


1687 


1901 


Kr,m* (MeV) 


832 


1073 


C(7tk;°)i,m* 


4.46 


2.84 




-3.41 


-0.13 


C (nR\ l ) ,M* 


16.8 


13.3 


C {-xRf) 2> M* 




0.15 




-0.76 


-10.0 


C^Rf^M* 




-0.17 


C(7Tiif )l,Af* 


10.4 


7.37 


C (ttR1 ) 3 ,M* 




-0.06 



TABLE V. Masses (M^,), cutoffs (A^jvf*), and couplings {Cr^ R Li\ M ») of the bare M* = 
7T2(1670). For the description of the table, see the caption of Table. HVl 

vr 2 (1670) 



Unitary model Isobar-fit model 



M°p (MeV) 
Kr,m* (MeV) 

C(7rR0°) 2 ,M* 

C(irR° )2,M* 
C( 7T Rl 1 ) 1 ,M* 


1877 

874 
0.67 
0.99 

-2.21 


1912 
885 
0.55 
0.97 
-1.67 


C {-kR\ 1 )z,M* 
^{nR^)i,M* 


0.50 


3.58 


C (-KRf) s ,M* 
C (-kR\°) ,M* 


-12.2 


-11.3 


C(7Tiif ) 2 ,M* 






CfrRf^M* 







B. Z-diagram effects on the pole positions of ai(1260), 7^(1670), and 7^(2100) 

In Sec. UVA2| we solved Eqs. §U]) and to fit the pole positions for ai(1260), 7r 2 (1670) 
and 7T2(2100) listed in PDG. Our fitted values are shown in the row labeled as "With Z" of 
Table IVIIi When the Z-diagram mechanisms are turned off, which is achieved by replacing 
the dressed vertex function F^ cR \ l>M * with the bare Ft cR \ l>M * in calculating T,m* (E) [Eq. fl51~j) ]. 
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TABLE VI. Masses (Mjjp), cutoffs (AkR,m*)i an d couplings (C^ r li^^ m , 
7r2(2100). For the description of the table, see the caption of Table. HVl 



of the bare M" 



vr 2 (2100) 



Unitary model 



Isobar-fit model 



M° M , (MeV) 
Kr,m* (MeV) 



C, 



'(ttR™ 
frR? 



c, 



2, M* 
2> M* 

i,M* 

3, M* 

1, M* 

3, M* 

a,M* 

2, M* 

4, M* 



2189 
876 
-0.93 
-0.001 
2.45 

0.51 

-11.9 



2280 
970 
-0.70 
-0.34 
1.65 

0.41 

-10.1 



TABLE VII. Pole masses of ai(1260), vr 2 (1670), and vr 2 (2100). Here, "With Z" denotes the results 
of the full unitary model, while "Without Z v denotes the results in which the Z diagrams are 
turned off from the full unitary model. 







Pole masses (MeV) 






ax (1260) 


tt 2 (1670) 


7T 2 (2100) 


With Z 


1230 - 213 i 


1672 - 130 i 


2090 - 313 i 


Without Z 


1122 - 148 i 


1661 - 127 i 


2044 - 398 i 



the solution of Eq. ( 1501) becomes the values shown in the row labeled as "Without Z" of 
Table IVHl 

Comparing the two rows in Table IVHl we see that the Z-diagram mechanisms can change 
the pole positions significantly. In particular, the imaginary parts can be changed by 65 MeV 
for ai(1260) and 85 MeV for 7r 2 (2100) Accordingly, we expect that the extracted residues will 
also be significantly changed. The extraction of the residues for unstable particle channels 



is nontrivial, and is still being investigated, as explained in Ref. [48]. We thus do not have 



results for the Z-diagram effects on the branching ratios in this work. 



C. Z-diagram effects on Dalitz plots 

1. D° vr+TT^vr 

As discussed in Sec. IIV A 21 we only include the bare D° — > irp vertex with ( J p = 0~, T 



0, / = 1) in this calculation, as guided by the analysis by the BABAR Collaboration [19 
the 7r/o and tt/2 channels are coupled only through the final state interaction. The Dalitz 
plot calculated from our unitary amplitude T 7r + 7r - 7r o j £)o(£' = M^o = 1865 MeV) from using 
Eqs. f lT9|) and (j4~3]) is shown in Fig. H] (left panel). With an overall normalization factor, 
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FIG. 4. (Color online) (left) Dalitz plot of D° 
distributions with and without the Z graphs. 



7r + 7r 7T° decay; (right) Ratio of Dalitz plot 



the pattern of our Dalitz plot is similar to BABAR's data 19j. The sharp peaks (darker 
red) near the edges of distributions are due to the formation of a p resonance during the 
3-7r propagation. The almost empty center part is due to the destructive interferences along 
the symmetry axes, supporting the assumption that the T = Tip channel dominates the 
decay (r32j . 

With the same parameters and overall normalization factor, we then calculate the Dalitz 
plot from T^kar^ D0 (£■) us ing Eqs. ( jl~9l) and ( 1491) which does not include Z-diagram mech- 
anisms, namely, keep only the first term in Fig. [2J In the right panel of Fig. HI we show 
the ratios between the results obtained from calculations with and without the Z diagram. 
Clearly, the Z-diagram mechanisms considerably change both the magnitudes and the shape 
of the Dalitz plot. In most of the area, the ratios (measured by the darkness as indicated on 
the right y axis of the figure) are about 1.6. To see this more clearly, we show in Fig. [5] the 
double differential decay width distribution, d 2 Y / (dM% +n0 dM%- w0 ) defined in Eq. ( 1B2j) . at 
M^ +7r o = 0.3 GeV 2 . We see that at the p resonance peaks, the magnitudes can be enhanced 
by a factor of about 1.5 when Z-diagram mechanisms are included to satisfy the three-body 
unitarity. 

Our results shown in Figs. H] and [5] indicate the need to reanalyze the D-meson decays 
data, with the three-body unitarity taken into account, to assess the results, such as CKM 
matrix elements, obtained with the isobar-model analysis [l8-23 . 



2. ai(1260),7r 2 (1670),7r 2 (2100) 



7T 7T 7T 



The decays of these three mesons have been analyzed by using the isobar models. Our 
objective here is twofold. First we want to examine the Z-diagram effects on the Dalitz 
plots. Second, we regard the Dalitz plot generated from our unitary model [Eq. ( I4"3"j) ] as the 
data, and fit them with the isobar model [Eq. (|49p ]. We refer to it as the isobar-fit model. In 
this way, we have the two models that reproduce the same Dalitz plot. However, the decay 
amplitudes from the two models are not necessarily the same, which we will examine. This 
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FIG. 5. (Color online) Double differential decay width distribution [Eq. (|B2|) ] of D° — > tt + tt~tt° 
decay at = 0.3 (GeV) 2 . The red solid curve is from the full unitary model while the blue 

dashed curve is obtained by just turning off the Z diagrams in the full model. 
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FIG. 6. (Color online) (left) Dalitz plot of oi(1260) 
distributions with and without the Z graphs. The unit is GeV 
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examination is particularly interesting in the context of the extraction of the CKM phase 7 
from B and/or D decays. This is because the extracted 7 depends on the decay amplitudes, 
particularly on its phased. Thus the difference in the decay amplitude between our unitary 
model and the isobar-fit model does matter. We examine this for the strong decays of a± and 
7T2, which is suggestive enough for the extraction of 7 with the isobar model analysis. Also, 
we study how well the three-body unitarity is satisfied in the isobar-fit model. Our unitary 
model satisfies it by definition (explicitly shown numerically later). A large violation of the 
unitarity raises a concern about the reliability in extracted quantities with the isobar-model 
analysis. 
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3 There exists an alternative approach in which 7 can be determined model-independently 
from data, provided a large dataset is available. A feasibility study jo^ showed that, with a dataset 
available in the near future, the precision of 7 extracted with this approach is comparable to that obtained 
with the isobar-model analysis. Future high statistic experiments (super B factory, LHCb) make this 
approach very interesting. 
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FIG. 7. (Color online) (left) Dalitz plot of vr 2 (1670) ->■ vr+vr-vr 
distributions with and without the Z graphs. The unit is GeV" 3 . 
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FIG. 8. (Color online) (left) Dalitz plot of vr 2 (2100) ->■ 7r + 7r~vr°; (right) ratio of Dalitz plot 
distributions with and without the Z graphs. The unit is GeV~ 3 . 



We calculate the Dalitz plots at E = Re[£? po i e ], for which we use -E po i c of the full model 
listed in Tabic IVIIl The results from our full model [using Eq. (l43j) ] are shown in the left 
panels of Fig. O for ^(1260), Fig. for tt 2 (1670), and Fig. for vr 2 (2100). We see that they 
have rather complex structure. This is due to the resonances of tttt scattering implemented 
in the tt-R Green function [Eq. (j41j) ]. and also interference among them as a consequence 
of summing coherently irR — > mrn partial- wave amplitudes as calculated in Eq. (143j) . For 
example, two bands on M^ +7r0 ~ 0.6 GeV 2 and M^ n0 ~ 0.6 GeV 2 in Figs. are due to 
the p(770) resonance in the p-wave tttt scattering. A gap in M 2 +7r0 + M^-^o ~ 0.6 GeV 2 
in Fig. [6] is due to the /o(980) resonance and opening of the KK channel. The Z-diagram 
effects are rather different in different parts of the Dalitz plots. This can be seen from the 
ratios between the Dalitz plots calculated with [Eq. (j4"3]) ] and without [Eq. (l4"9|) ] Z diagram, 
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TABLE VIII. The pole masses and total widths decaying to 7r 7r + 7r° states (r^+^-^o) of ai(1260), 
vr 2 (1670), and vr 2 (2100). "Full" ("Isobar fit") is the results of the full unitary (isobar-fit) model. 



M* 


ai(1260) 




vr 2 (1670) 




7T 2 (2100) 






Pole mass 


n 


Pole mass " 


L 7T+7T - 7T° 


Pole mass ] 


L 7T+7T - 7T° 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


Full 


1230 - 213z 


375.4 


1672 - 130i 


157.8 


2090 - 313i 


219.0 


Isobar fit 


1230 - lOOi 


371.9 


1672 - 97i 


151.2 


2090 - 261* 


217.0 



as shown in the right-hand sides of Figs. | 

To see the Z-diagram effects more clearly, we show in Fig. [5] the double differential decay 
width distributions [Eq. HB2j) ] for the decays of these three mesons at typical kinematics. By 
comparing the red solid curves and blue dashed curves, we see that the Z-diagram effects 
can significantly reduce the cross sections; in particular in the regions near the resonance 
peaks of irir scattering. 

Now we examine differences between the unitary and isobar models if both fit the same 
Dalitz plot. We treat the Dalitz plots in the left sides of Fig. |6]-Fig. [8] as the data in fits 
using the isobar model [Eq. (}4l?]) ] by adjusting all the available coupling constants, cutoffs 
of the vertex Fr cR \ uM * of Eq. ( H7|) . In the fits, we assign either 5% error for each point of the 
Dalitz plot larger than 0.005 GeV -3 , or error of 0.005 GeV -3 otherwise. We are able to get 
reasonably good The resulting parameters for the isobar-fit model are rather different 
from the unitary model as shown in Tables IIV1IVH The quality of our fits can be seen by 
comparing the dotted curves and the solid curves in Fig. [9j Accordingly, the decay widths 
to 7T + 7r~7T channel calculated from two models using Eq. (I31~j) (keeping only abc = 7T + tt^7t ) 
agree well, as seen in the third, fifth and seventh columns of Table IVIIII However, we see in 
the second, fourth and sixth columns of Table IVIIII that the resonance pole positions from 
resulting isobar-fit models differ significantly from those of the unitary model from which 
the Dalitz plot data are generated. Their imaginary parts can differ by more than 100 MeV 
for ai(1260) and about 50 MeV for 7r 2 (1670) and 7r 2 (2100), indicating a large violation of 
the three-body unitarity. Note that the bare M * mass [Mj^* in Eq. (126]) ] does not enter the 
calculation of the Dalitz plots. Thus we choose M° 7 * for the isobar-fit model so that the 
real part of the pole is the same as that of the unitary model. 

Some hadron models predict that the hybrid mesons can have quite different branching 
ratios from those of the ordinary mesons with radial excitations of the quark-antiquark 



pair 66]. Thus, the parameters, C^r^m*-, m Tables ITVWll can provide important information 
to distinguish the hybrid and/or exotic mesons from the ordinary mesons. The significant 
difference in the parameters between the unitary and isobar models indicates that we should 
use a unitary model to analyze the Dalitz plot distributions to search for the exotic mesons. 

As discussed above, the importance of using an unitary model can be seen more clearly 
in comparing the M* — > irR amplitudes predicted by the two models. This is shown in 
Fig. [10J The M* — > ttR amplitudes generated from unitary model must be complex because 
of multiple scattering due to Z-diagram mechanisms, while those from the isobar model 

4 In most Dalitz plot analyses with the isobar model, a (constant) nonresonant background amplitude with 
adjustable strength is included. Also, the M* — > ttR couplings arc in general complex in the isobar model 
in order to partially take account of the missing interaction between the spectator and the paired mesons. 
Because we obtained fits good enough for the following discussion, we do not include these degrees of 
freedom. 



22 




FIG. 9. (Color online) Double differential decay width distributions [Eq. (|B2p ] of ai(1260) — > 
tt+tt-tt (top), of tt 2 (1670) tt + tt~tt° (middle), of vr 2 (2100) -> TT+vr-vr (bottom). The red solid 
curves are from the full unitary model while the blue dashed curves are obtained by just turning 
off the Z-diagrams in the full model. The black dotted curves are from the isobar-fit model (see 
the text for definition.). 
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FIG. 10. (Color online) The partial-wave M* — > cR vertices. The dressed (bare) vertices of the 
full (isobar-fit) model defined in Eq. ([36]) [Eq. (fTF|) ] are shown by the red solid (blue dashed) 
curves. The upper (lower) figure shows the absolute value (phase) of the vertices in arbitrary 
scale (degrees). The figures on the left, middle and right columns are for ai(1260), 7^(1670) and 
^(2100), respectively. The label in each figure specifies the cR state. We show the result for the 
bare R = /o or p of the lowest bare mass and c = tt, and the relative orbital angular momentum 
of the cR is specified by Ir c . 



can be chosen to be real (cf. footnote H]). Their differences in real parts can also be very 
different in some regions. The difference in the phase is more apparent. Even though 
nonzero phase could have been used in the isobar-fit model, as has been done in most 
isobar-model analyses, the rather large dependence on the kinematics, which reflects the 
three-body unitarity, is beyond the capability of the isobar model to simulate. As we have 
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TABLE IX. Comparison between the total decay width of bare M* (X^-k+ttKr) and twice of the 
imaginary part of the M* self-energy (— 2Im[£jw*]). Both of r 37r+7rA -^ and Em* are calculated at 
E = Re[i2p i e ], where E po \ c is a pole mass of a physical M* listed in Table IVIII1 "Full" ("Isobar 
fit") is the results of the full unitary (isobar fit) model. 



M* 


ai(1260) 


vr 2 (1670) 


vr 2 (2100) 






— 2Im[£jy/*] 


^3n+nKK 


-2Im[£ M *] 




-2Im[E M .] 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


Full 


379.8 


379.7 


234.7 


234.8 


413.7 


413.4 


Isobar fit 


378.4 


266.2 


227.4 


198.7 


434.5 


379.1 



noted, the phases of these amplitudes are crucial in using .D-meson decays to determine the 
phase 7 of CKM matrix elements as a way to find physics beyond the standard model. The 
previously extracted 7 from B T — > D° (or D°) K T — > (Kgii + 7i~)K T has the uncertainty 



from the isobar model fitted to the D-decays. It is estimated to be 8.9° for Belle [23], 3° for 



BABAR j2l|. Considering the difference in the phase, typically of 10° ~ 20° level, between 
the unitary and isobar models, it would be highly desirable to analyze the data with the 
unitary model. 

Finally, let us examine the extent to which the three-body unitarity is satisfied by each 
model. We can examine this using Eqs. fj3~T|) and (132]) which are satisfied by a unitary model. 
Within our current model developed for ai(1260), 7r 2 (1670) and 7r 2 (2100) decays, the total 
decay width {T 37T+7rK ^) is the sum of M* — > 3n (T 37T ) and M* — >• txKK (T wK k) widths. 
In the row labeled by "Full" of Table IIX[ we can see that the unitarity relation is satisfied 
within the numerical precision, as it should be. On the other hand, for the isobar-fit model, 
the unitarity is rather badly violated as seen in the fifth row of Table IIXI which raises a 
concern about the reliability of results obtained with the isobar model. 



V. SUMMARY AND OUTLOOK 



Starting with a model Hamiltonian with vertex interactions f a b,R and T c r m* and two- 
body interactions f c '_R', c i?, where R and M* are the bare one-particle states and a, b, c are 
light pseudoscalar mesons (ir, K, etc.), we have developed a unitary coupled-channels model 
for three-mesons decays of heavy mesons and excited meson states. By fitting the empirical 
amplitudes for meson-meson scattering such as nn tttt, the vertex interactions f a b,R, which 
can generate resonances R in meson-meson scattering, are determined and used to predict 
the one-particle-exchange Z-diagram mechanisms Z c iRi fi n(E). The scattering amplitudes 
T c >r', c r(E) are then calculated with Z c iRr tC n(E) by solving a set of coupled-channels equations 
with the three-body unitarity condition satisfied exactly. The final state interactions of 
three-mesons from the decays of heavy mesons are then calculated from T c / # iC r(E) . In the 
absence of the Z-diagram mechanisms, our decay amplitude is reduced to a form similar 
to that used in the isobar-model. This allows us to investigate the extent to which the 
commonly used isobar-model analysis is valid in extracting the properties of heavy mesons 
from the Dalitz plots of the measured three-mesons distributions. For strong decays of a 
heavy meson M*, we present formula and procedures for investigating the importance of 
three-meson interactions in determining the resonance pole positions on complex-energy 
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Riemann surface. 

The model has been applied to investigate three-pions decays of ai(1260), 7r 2 (1670) and 
7T2(2100), and D° mesons. It was found that the Z-diagram mechanisms can change signif- 
icantly the magnitudes and shapes of the Dalitz plots. For D° — > 7r + 7r~7r°, the changes in 
magnitudes can be a factor of about 1.6 in most of the phase space. For ai(1260), 7r 2 (1670) 
and 7^(2100), the changes are about a factor of 1.3 ~ 1.6 in magnitudes in the regions where 
meson-meson resonances /o(600), p(770) and /2(1270) dominate. We have also examined 
differences between the unitary and isobar models, both of them produce the same Dalitz 
plot. We have demonstrated that decay amplitudes from the two models are significantly 
different, particularly in the phase. A proper estimate of the phase is particularly important 
for extracting the CKM phase 7 from data of B T -> D° (or D°) K T -> (K%7i + n-)K T 
for which the D decay Dalitz plot is analyzed with a model. We have also shown that 
the three-body unitarity is rather largely violated in the isobar model. Finally, the result- 
ing bare parameters which can be interpreted as characterizing the "intrinsic" quark-gluon 
substructure of heavy mesons are also very different between the unitary and isobar models. 

Our results strongly indicate the need for reanalysis of the D-meson decays using a uni- 
tary model to assess the results, such as the CKM matrix elements, obtained with the 
isobar- model analyses [l8-23|]. It is also important to reanalyze the three- meson decays of 



all heavy mesons listed by PDG as a necessary step for establishing meson spectroscopy 
and exploring the hybrid or exotic mesons in the near future at JLab, GSI, and other pos- 
sible facilities. While the model presented in this work is more advanced than the models 
used in the previous analyses of three-meson decays processes, improvements are needed to 
make quantitative progress. We need to include the data associated with KK channels in 
determining our parameters. The approach for extending our formulation to include effects 
due to four-pions channels, which are considered to be important for determining scalar- 
isoscalar (L = I = 0) resonances, should be developed. For 5-meson decays, an appropriate 
theoretical approach must be developed to describe irir amplitudes at high energies where 
no data is available. Finally, our formulation is derived from applying a unitary transfor- 
mation 67, 68| to a Hamiltonian defined within the relativistic quantum field theory. As 
discussed in Ref . 39| , this method, as well as many well-studied three-dimensional reduction 



methods [69|, is needed to derive tractable reaction models for solving complex reactions 
involving many channels and three-particle final states, with the unitarity maintained. Nev- 
ertheless, accuracy of these approximations should be investigated in the future. 
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Appendix A: Relation between E^* and M* decay amplitude 



Here we derive Eq. (|32|) . In the course of the derivation, we will see that Eq. (132]) holds 
true only when the T matrix (X") satisfies the unitarity relation. We start with the "decay" 
width T t A 4* (E) of the bare M* (not the physical resonance state) defined by 



t * (E) =2nJ2 KE -E a -E b - E c ) \T abc , M * (£) | 5 



abc 



The amplitude T abCjM * is defined in Eq. (TT3|) . and it can be written as 

T abc , M * = (abc\(l + T'G )H'\M*) , 



(Al) 



(A2) 



where Go is the free Green function, and the reaction T-matrix (T") has been defined in 
Eq. dSJ). By using Eqs. (lAlj) and (IA2[) together with the unitarity relation 



T - T /f = -27riT'5(E - H Q )T'\ 

and the equality 

G -G = -2-Ki5(E-H ) : 

we arrive at Eq. (1521) as 

r t A £ (E) = 2tt £ S{E -E a -E b - E c ) \T abcM , (E)\ 2 

abc 

= 2tt]T<5(£ -E a -E b - E c ) (M*\H'(1 + G T't)|a6c) 



(A3) 



(A4) 



(the 



= 2ir(M*\H'(l + GlT'i)5(E - H )(l + T'G )H'\M*) 

= i [(M*\H'G Q (1 + T'G Q )H'\M*) - (M*\H'{1 + G\T^)GIh'\M*) 

= -2 Im [(M*\H'G (1 + T'G )H'\M*)} 

= -2 Im [E M . (£)]• 

In the last step we have used the definition of E m* (-E 1 ) given in Eq. ( 1271 ) . 



(A5) 



Appendix B: Dalitz plot 

Here we summarize the formulae for calculating Dalitz plots. The differential decay width 
of a heavy meson M* (E: mass; Sm*' spin) at rest decaying to three (pseudo) scalar mesons, 
M*(0, E) — > a(p a ) + b(pb) + c(p c ), where we consider only (pseudo) scalar mesons for the 
final state, can be expressed as 

, r 1 d 3 Pa d% a*Pc B ^ . , 2 

M * 2E {2nf2E a {p a ) (2ir)*2E b (p b ) (2tt) 3 2£ c ( Pc ) 2S M * + abcM * 1 

x (2vr) 4 5(£ - £ a (p a ) - £ 6 (p&) - E c {p c ))5 3 (0 -Pa-Pb- Pc), (Bl) 

where M. a bc,M* is the invariant amplitude of the decay; B is the Bose factor for the final 
mesons. For example, when the three final mesons are identical, B = 1/(3!). With a 
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variable transformation, we obtain the double differential decay width distribution (Dalitz 
plot density) for the unpolarized decay given by 

dm 2 ab dm 2 c = (27r)3 32£3 2S M * + 1 £ \ M ^ M '^ ( B2 ) 

M* 

where m a b (m bc ) is the invariant mass of the ah (be) pair. The invariant amplitude is related 
to the decay amplitude defined in Eq. ffT9l) and Eq. (123]) or Eq. (113"]) by 



MabcM* = -(2n)W2Eyj2E a ( Pa )yj2E b (p b ),j2E e (p B )T abe ,M*. (B3) 

The meson labels a, 6, c specify the momentum (p x ), the mass (m x ), and the isospin (t x ) of 
the meson x = a,b, c. 

Next, we summarize relations between kinematic variables. For a given value of m 2 b , the 
range of m\ c is determined by its values when p b is parallel or antiparallel to p c : 



(Omax = {El + E* c f - JEf -m 2 - JEf - m 2 , 



lUn = (E* b + E* c f - JEf -m 2 + JEf -ml) , (B4) 



with 

b = 2m~^ mab ~ ma + mh >' 

E * = ^f h ( E2 - m2 a b + m c)^ ( B5 ) 

being the energies of the particles b and c in the center-of-mass frame of the ab pair, respec- 
tively. For a given set of m a b and m& c , the momenta of the final particles are 







Pa = 


2E^ 






Pc = 


2E^ 


Pb = 


A 2 - 




i 


®ab = 


2p a Pb 



'b- 



(E-E a -E b f-m 2 c -p 2 a -p 2 b \, (B6) 

where 6 a b is the angle between p a and p b . Taking p a on the xz-plane, we have 
Pa = Pa(sm9 a ,0,cos9 a ), 

p b = Pb(cos 6 a sin 9 ab cos 4> ab + sin 6> a cos 9 ab , sin 6> ab sin <f) ab , 

- sin 6> a sin 9 ab cos </> ab + cos 9 a cos ab ) , 
Pc = ~Pa - Pb, (B7) 

where <j) a b is the azimuthal angle of p b > which is obtained by rotating p b around |/-axis by 
— 9 a . To calculate the differential decay width for the unpolarized decay [Eq. (IB2I) ]. one may 
set 9 a = and (p ab = 0. 
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Appendix C: Z-diagrams 



1. Definition 



The matrix element of the Z-diagram for a transition process, R(—p c )+c(p c ) —> R'(—p d ) + 
c'(p d ), is given by 

{R'(-p d , s z R ,,t z R ,); c'(p d , 0, t z c ,) \Z C "(E) \R(-p c , s z R , f R ); c(p c , 0, t z c )) = 

(R'(-Pc',s z RI ,t z R ,)\f R/ y /c \c(p CJ 0, **); c"(p d >, 0, t z „)} 

1 

V 

E - E c (p c ) - E d (p d ) - Ec»(p c n) + ie 
x (c' (p d , 0, #); c" (p c «, 0, ^0l/cc",fl|i2 (-Pc, 4, **)>• (CI) 

Here c" is the exchanged meson; s R (t z R ) is the ^-component of the spin (isospin) of the 
particle R; t z is the ^-component of the isospin of the particle c; p c « = — p c — p d . The 
vertices are expressed by 



(c\p d , 0, t z d ); c"(p d ,, 0, t z c „)\f dd ,, R \R (-p c , s z R , t z R )) = 

Jr(Pc',Pc", qc)(t c 't z d t c /d z d ,\t R t z R )Y SRS z R (q c )f c i c // >R (q c ), 
(R' (-p c r, s R ,,t z R ,)\f R , tC „ c \c {p c , 0, t z )\ c" (p c „, 0, £,)) = 

JR' (Pc, Pc", Qc'){tc"t Z c ^t c t Z c \t Rl t Z R ,)Y* RiS ^ (q c ')fR',c"c(q C ' 



(C2) 



(C3) 



The above equations are the same as Eq. (JH]), and / is related to the ttti model [Eq. fl35l) ] 
through Eq. fl39l) . Here q c is the meson momentum in the center of mass of the two-meson 
(c'c") subsystem from the R(—p c ) —> c'{p d ) + c"(p c ») decay The factor J R appears as a 
result of the Lorentz transformation of the vertex, and is given by 



jR(Px,Pc",q v ) 



E x (q y )E dl (q y )m R 



^ E x (p x )E d/ (p c n)E R {py) 
with x, y = c or c'. Using the Lorentz transformation, we have 

Qc = Pc' ~ p(-Pc,Pc')Pc = K c p c + \ c p d , 
-q c > =Pc~ p(-'Pc',Pc)Pc' = K d p c + X d p d 



where 



P{P,Pz 



and 



p-p. 



£(P, p x ) + E x {p x ) + E d , {\P-p x \) 



- E>x{Pi 



E x {p x ) + E d ,{\P - p x \) 



P 2 . 



(C4) 



(C5) 
(C6) 



(C7) 



(C8) 



The signs attached to p c and p d in Eqs. (1C5I) and (1C6I) . and the ordering of isospins in the 
Clebsch-Gordan coefficients in Eqs. (IC2I) and (1C3I) matter for the phases of the p -H- titi 
interaction, and are taken consistently with the pmr interaction Lagrangian. To take these 
phases appropriately is important for giving the correct phases to the amplitudes. 



29 



2. Partial-wave decomposition of Z potential 

We define the partial-wave expansion of Eq. ( 1C1I) as 

(R'(-p c/ , s R „t RI ); c'(p c ,, 0, t z ,)\Z c "(E)\R(-p c , s z R , t z R ); c, (p c , 0, t z c )) = 

EE E (t R t R u c \TT*)(t R d R d^\TT*) 

TT Z JJ Z VV z ll z 

x(ll z s R s z R \JJ z }{l'l /z s R ,s z R ,\JJ z } 

x Y VjV . (pcf)Y^ lt (Pc)Zp R ^ ){cR)l (p c i,p c 

Performing some manipulations, we obtain 

(-!)*="-*«+*=' J(2t R + l)(2t R , + l)W(t c tjy^; t c »T) 



(C9) 



yc",JPT 



x(-l) s «'V(2 Si? / + l)(2s fl + 1)(2/' + 1)(2Z + 1) 



X (2j + l)(2L' + l)(2L" + l 

l a ,l b ,L',L",j 



(25^ + 1)1(23^+1)! 



X 



\(2l a )\(2s R ,-2l a )\(2l b )\(2s R -2l b )\ 
V l a L'\( I l b L" \(s R - l b j L'\( s R , - l a j L" \ 



/ \ 











x j /' i. L' If / Z 6 L" 1 f L" - l h J 1 FWb 



(C10) 



where we have used (— iy+ l '+ s R+ s R' = 1 from the parity conservation. We have introduced 
j 



F) a ' h and B defined by 



1 r 1 



dx 



BP 3 (x) 



3 2 7-i E — E c (p c ) - E c , {p d ) - E d , {p c „) + le ' 

B = J{p c , p c " , q c > ) f S J*R> (q c ')J (Pc> , Pc" , g c ) (?c) 

x(A c ^ c /Q (^ c ) s «'- /o (A c p c s «" Zb (^ c )' 6 (^)" Sii (^)"^' ; 



(Cll) 



(C12) 



where Pj (x) is the Legendre function of the degree j . 
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